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Abstract. The Smoluchowski equation is a system of partial differential equations modelling the 
diffusion and binary coagulation of a large collection of tiny particles. The mass parameter may be 
indexed either by positive integers, or by positive reals, these corresponding to the discrete or the 
continuous form of the equations. In dimension d > 3, we derive the continuous Smoluchowski PDE 
as a kinetic limit of a microscopic model of Brownian particles liable to coalesce, using a similar 
method to that used to derive the discrete form of the equations in [4J . The principal innovation is 
a correlation- type bound on particle locations that permits the derivation in the continuous context 
while simplifying the arguments of |lj . We also comment on the scaling satisfied by the continuous 
Smoluchowski PDE, and its potential implications for blow-up of solutions of the equations. 



1 Introduction 

It is a common practice in statistical mechanics to formulate a microscopic model with simple 
dynamical rules in order to study a phenomenon of interest. In a colloid, a population of compar- 
atively massive particles is agitated by the bombardment of much smaller particles in the ambient 
environment: the motion of the colloidal particles may then be modelled by Brownian motion. 
Smoluchowski's equation provides a macroscopic description for the evolution of the cluster densi- 
ties in a colloid whose particles are prone to binary coagulation. Smoluchowski's equation comes in 
two flavours: discrete and continuous. In the discrete version, the cluster mass may take values in 
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the set of positive integers, whereas, in the continuous version, the cluster mass take values in R^. 
Writing /„(x,t) for the density of clusters (or particles) of size n, this density evolves according to 

(1.1) ^ = d{n)AUx,t)+Ql{f){x,t) - Qlif){x,t), 
where 

/•n 

(1.2) Ql{f)= P{m,n-m)fmfn-mdm, 

Jo 

and 

(1.3) Q^_{f)=2 P{m,n)fmfndm, 

Jo 

in the case of the continuous Smoluchowski equation. In the discrete case, the integrations in (|1.2p 
and (|1.3|) are replaced with summations. 

In m and [5], we derived the discrete Smoluchowski equation as a many particle limit of a 
microscopic model of coagulating Brownian particles. (See also [7], [9] and [2] for similar results.) 
The main purpose of the present article is the derivation of (1.1) in the continuous case. We 
introduce a simpler approach to that used in [3] and [5] . We will present a robust argument that 
allows us to circumvent some induction-based steps of [1] and [3] (which anyway could not be 
applied in the continuous case). As such, an auxiliary purpose of this article is to present a shorter 
proof of the kinetic limit derivations of Smoluchowski's equation given in [3] and [5]. The main 
technical tool is a cor relation- type bound on the particle distribution that seems to be applicable 
to general systems of Brownian particles. To explain this further, we need to sketch the derivation 
of Smoluchowski's equation and explain the essential role of the correlation bounds. 

The microscopic model we study in this article consists of a large number of particles which 
move according to independent Brownian motions whose diffusion rates 2d{m) depend on their 
mass m E (0, oo). Any pair of particles that approach to within a certain range of interaction are 
liable to coagulate, at which time, they disappear from the system, to be replaced by a particle 
whose mass is equal to the sum of the masses of the colliding particles, and whose location is a 
specific point in the vicinity of the location of the coagulation. This range of interaction is taken 
to be equal to a parameter e, whose dependence on the mean initial total number N of particles is 
given hy N = k^Z for a constant Z, where 



e^-'^ if d > 3, 
|loge| ifd = 2. 



This choice will ensure that a particle experiences an expected number of coagulations in a given 
unit of time that remains bounded away from zero and infinity as is taken to be high. 

Our main result is conveniently expressed in terms of empirical measures on the locations Xi{t) 
and the masses mi{t) of particles. We write g{dx, dn, t) for the measure on M*^ x [0, oo) given by 

g{dx,dn,t) = k^'^^6(^^(t),m,(t)){dx,dn). 

i 
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Our goal is to show that, in the low e limit, the measure g converges to fn{x,t)dx dn, where /„ 
solves the system (jl.ip . The main step in the proof requires the replacement of the microscopic 
coagulation propensity a{n, m) (that we will shortly describe precisely) of particles of masses n and 
m with its macroscopic analogue (3{n,m). The main technical tool for this is a correlation bound 
which reads as follows, in the case that the coefficient d{m) is non-increasing in m: 

/CO d 
K{xi,{t),..,Xi,{t))Y[d{mi^{t))^mi^{t)dt 

ii,...,ik i'=l 

k 

< const. E ^ k{xi,{0),mi^{0),...,Xi^iO),mi,{0))Yld{m,M)^miAO)- 

h,...,ik r=l 

Here, E denotes the expectation with respect to the underlying randomness, K : (M*^)^ ^ M is 
any non-negative bounded continuous function, and K = — ^^(mjJA^:^^ + • • • + (i(mij.)A2,.-^^ K. 
We refer to Section 4 for the corresponding correlation inequality when the function d(-) is not 
non-increasing. 

In fact, we need ()1.4p only for certain examples of K with k = 2,3 and 4. It was these examples 
that were treated in and [5] with rather ad-hoc arguments based on an inductive procedure on 
the mass of the particles. Those arguments seem to be specific to the discrete case and cannot 
be generalized to the continuous setting. Moreover, the bound (II. 4p implies that the macroscopic 
particle densities belong to for given p >2, provided that a similar bound is valid initially. This 
rather straightforward consequence of (jl.4p is crucial for the derivation of the macroscopic equation. 
The corresponding step in [1] and [5] is also carried out with a method that is very specific to the 
discrete case and does not apply to the continuous setting. This important consequence of ()1.4p 
simplifies the proof drastically and renders the whole of section 4 of [1] redundant. 

We state and prove our results when the dimension is at least three. However, our proof for the 
correlation bound ()1.4p works in any dimension, and an interested reader may readily check that, 
as in this article, the approach of [5j may be modified to establish Theorem 1.1 in dimension two. 

We continue with the description of the microscopic model and the statement of the main result. 

As a matter of convenience, we introduce two different microscopic models, that differ only in 
whether the number of particles is initially deterministic or random. We will refer to the model as 
deterministic or random accordingly. In either case, we define a sequence of microscopic models, 
indexed by a postive integer N. 

A countable set / of symbols is provided. A configuration q is an M.'^ x (0, oo)-valued function 
on a finite subset Jq of /. For any i € /q, the component q{i) may be written as {xi,mi). The 
particle labelled by i has mass rrii and location Xj. 

In the deterministic case, the index of the model specifies the total number of particles present 
at time zero. Their placement is given as follows. There is a given function h:W^x (0, oo) [0, oo), 
with hn{x) := h{x,n), where f^a h{x,n)dxdn < oo. We set Z = f^dhn(x) dxdn e (0, oo) 
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and choose points in (0, oo) x independently according to a law whose density at (x, n) is 
equal to hn(x)/Z. Selecting arbitrarily a set of N symbols {ij : j € {1, . . . ,N}} from I, we define 
the initial configuration q(0) by insisting that qi- (0) is equal to the j-th of the randomly chosen 
members of (0, oo) x W^. 

In the random case, the index gives the mean number of initial particles. We suppose given 
some measure 7Ar on positive integers that satisfies K{y]y) = N and Var(7jv) = o{N'^). The initial 
particle number, written M, is a sample of 7iv. The particles present at time zero are scattered in 
the same way as they are in the deterministic case. The subsequent evolution, whose randomness 
is independent of the sampling of Af, is also the same as in the deterministic setting. 

To describe this dynamics, set a parameter e > according to = k^Z, as earlier described. 
Let F : {R.'^ x [0,oo)}-^ [0, oo) denote a smooth function, where its domain is given the product 
topology. The action on F of the infinitesimal generator L is given by 

(LF)(q) =AoF(q)+AeF(q), 

where the diffusion and collision operators are given by 

AoF(q) = d{m,)A^^F 

and 

(1.5) A,F{q) = 6-V(^i^)a(m,,,m,) 

l^^FiSl^d) + -^^FiSl^d) - F(q)" . 

Imi + nij '■' nii + rrij J 

Note that: 

• the function ^ : M'^ ^ [0, oo) is assumed to be Holder continuous, of compact support, and with 
/iRd Vix)dx = 1. 

• we denote by SjjCi that configuration formed from q by removing the indices i and j from Jq, 
and adding a new index from / to which SjjCi assigns the value {xi,mi + mj). The configuration 
Sfjq is defined in the same way, except that it assigns the value {xj,mi + rrij) to the new 
index. The specifics of the collision event then are that the new particle appears in one of the 
locations of the two particles being removed, with the choice being made randomly with weights 
proportional to the mass of the two colliding particles. 

Convention. Unless stated otherwise, we will adopt a notation whereby all the index labels 
appearing in sums should be taken to be distinct. 

We refer the reader to [1] and [10] for the reasons for choosing = e'^"^Z, the form of the 
collision term in (1.5), and the interpretations of the various terms. 
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Let us write Mzi^"^ x [0, oo)) for the space of non-negative measures tt on R'' x [0, oo) such 
that 

TT (r'^ X [0, oo)) < Z. 

This space is equipped with the topology of vague convergence which turns A4z into a compact 
metric space. We also write Al^^R'^ x [0,oo)^) for the space of non-negative measures n such 

that for every positive T, /U^M'^ x [0, oo) x [0, T]j < TZ, which is also compact with respect to the 
topology of vague convergence. This space has a closed subspace X which consists of measures 
such that x [0, oo) x [ti,t2]) < (^2 — ti)Z, for every ti < t2- As we will show in Lemma 6.2 

of Section 6, the space X consists of measures ii{dx,dn,dt) = g{dx,dn,t)dt with t i— > g{dx,dn,t) 
a Borel-measurable function from [0, oo) to Mz(R'^ x [0, oo)). We will denote by Pjv = the 
probability measure on functions from t G [0, oc) to the configurations determined by the process 
at time t. Its expectation will be denoted Ejv. Setting 

g%dx,dn,t) = e'^~'^^6(^^.^t),mi(t)){dx,dn), 

i 

the law of 

q H- > g'^{dx, dn, t)dt 

with respect to induces a probability measure on the space X. We note that, since the space 
X is & compact metric space, the sequence Ve is precompact with respect to the topology of weak 
convergence. 

For the main result of this article, we need the following assumptions on q;(-, •) and d{-): 
Hypothesis 1.1. 

• The diffusion coefficient d : (0, oo) (0, oo) is a bounded continuous function and there exists a 
uniformly positive continuous function : (0, oo) (0, oo) such that both (/){■) and (l){-)d{-) are 
non- increasing . 

• The function a : (0, oo) x (0, oo) — (0, oo) is a bounded symmetric continuous function satisfying 

a(n, m) 

sup sup 3 < oo, 

n<L m md{m)^ (p{m)^~^ 

for every L > 0. 
Remarks. 

• The condition that the function cf) : (0, oo) — > (0, oo) exist is rather mild and is satisfied if d{-) is 
non-increasing. This condition requires that heavier particles to diffuse slower which is natural 

from a physical point of view. In fact when d{-) is non-increasing, then we can simply choose 
(p{m) = 1. Also, if d{-) is non-decreasing, then the function (p exists and can be chosen to 
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be (j){m) = d{m)~^. From these two cases, we guess that the first condition is related to the 
variation of the function d(-). As we will show in Lemma 2.2 of Section 2, the existence of such 
a function (p is equivalent to assuming that the total negative variation of \ogd{-) over each 
interval [n, oo), n > 0, is finite. 

• We note that if the function d{-) is non-increasing, then the second condition for small m and 
n is equivalent to saying that a{m,n) < Cmin(m, n). However, when m and n are large, the 
second condition is satisfied if for example a(m, n) < Cmd{m)'^/'^nd{n)'^/'^ . In summary, the 
second condition is rather mild if m and n are large, but much more restrictive if both are small. 
Our stipulation that d be bounded is more restrictive in the case for values of its argument close 
to zero, since it is reasonable to assume that very light particles diffuse rapidly. 



We also need the following assumptions on the initial data h: 
Hypothesis 1.2. 



• Jp°° / hn{x)dxdn < oo. 

• * Xk ^ ^hcO^'^)^ fc = 2,3 and 4, where = nd{n)^~^ (p{n)^~^hn dn and Afc(x) = 

jKx)hiy)\x-yr^dxdy<oo 

where h = J^{n + l)hndn. 



Remark. Recall that if d{-) is non-increasing, then we may choose = 1. In this case, Hypothesis 
1.2 is satisfied if h e D L°°. 

To prepare for the statement of our main result, we now recall the weak formulation of the system 
(jl.ip . Firstly, recall that a non-negative measurable function / : M'^ x [0,oo) x [0,(X)) — > [0, oo) is 
a weak solution of (1.1) subject to the initial condition f{x,n,0) = hn{x), if for every smooth 
function J : M'^ x (0, oo) x [0, oo) — > M of compact support, 

f{x,n,t)J[x,n,t)dxdndt = / / hn{x)J[x,n,0)dxdn 

Jo J 

r* r Qj 
+ 11 I --—(x,n,s)f(x,n,s)dxdnds 

ft roc r 

+ / d{n)AJ{x,n,s)f{x,n,s)dxdnds 

Jo Jo Jr"^ 

ft fOQ fOQ f 

+ / / / / PiiT^, 'IT') f{x,n, s)f{x,m, s) 
Jo Jo Jo Jr'^ 

J{x, m, n, s)dxdndmds, 
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where 

J(x, m, n, s) = J{x, m + n, s) — J{x, m, s) — J(x, n, s). 
Following Norris [9] , we define an analagous measure- valued notion of weak solution. 

Definition 1.1 Let us write M[0, 00) for the space of non-negative measures on the interval [0, 00). 
We equip this space with the topology of vague convergenece. A measurable function / : M'^ x 
[0,00) M[0,oo) is called a measure-valued weak solution of (EiP ^Z; for each i > 0, the 

functions ge,he G Lj^^, where 

gi{x,t) = f{x,t,dn), hi{x,t) = / P{m,n)f{x,t,dn)f{x,t,dm), 
Jo Jo Jo 

and, secondly, 



00 



J{x,n,t)f{x,t,dn)dx = / / hn{x)J{x,n,0)dxdn 

rt rco r Qj 







yX,n,s)f{x,s,dn)dxds 

10 JO JR'^ c't 

(1.6) ~^ / d{n)AJ{x,n,s)f{x,s,dn)dxds 

Jo Jr'' Jo 

+ / / / / P{m,n)J{x,m,n,s)f{x,s,dn)f{x,s,dm)dxds. 
Jo JR'' Jo Jo 

Remark: The requirement ge,hi € Lj^^ is made in order to guarantee the existence of the integrals 
in (fLBD. 



We are now ready to state the main result of this article. 

Theorem 1.1 Consider the deterministic or random model in some dimension d > 3. Assume 
Hypotheses 1.1 and 1.2. If V is any limit point of T"^ , then V is concentrated on the space of 
measures g{dx, dn,t)dt = f{x,t,dn)dxdt which are absolutely continuous with respect to Lebesgue 
measure dx x dt, with f solving the system of partial differential equations (jl.ip in the sense of 
(1.6). The quantities (3 : (0, 00) x (0, 00) (0,oo) are specified by the formula 

f3{n,m) = a{n,m) / V{x) [1 -\- u{x;n,m)] dx, 



where, for each pair (n, m) G (0, 00) x (0, 00), u{-) = u{-; n, m) : M"^ — > (0, 00) is the unique solution 
of 



1.7 An X = ^. ; ' / . V{x) 

^ ' ^ ^ d{n) + d{m) ^ ' 

satisfying u(x) ^ as |xj 00. 



1 + n(x) 
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Remarks. 



• The continuity with respect to m and n and other important properties of u{-;n,m) will be 
stated in Lemma 4.2 of Section 4. In particular u G [—1,0], which implies that /? > because 
u is not identically zero. It follows from Lemma 4.2 that /3 is a continuous function. We also 
refer to the last section of [4j in which several properties of /? are established. In particular, it 
is shown that (3 < a and P{n,m) < Cap{K){d{n) + d{m)), where K denotes the support of the 
function V and Cap(i^) denotes the Newtonian capacity of the set K. (See [1] for the definition 
of Newtonian capacity.) 

• To simplify our presentation, we assume that all particles have the same "radius". However, 
in a more realistic model, we may replace £~'^V{e~^{xi — xj)) with e~'^V{e~^{xi — Xj);'mi,mj), 
where V{a]n,m) = (r(n) -\- r{m))~'^V{a/{r{n) +r(m))) and r(n) is interpreted as the radius of 
a particle of mass n. Our method of proof applies even when we allow such a radial dependence 
and we can prove Theorem 1.1 provided that r(n) = with x < {d — 2)~^ (when d > 3). In 
fact, we anticipate that, if x > (d — 2)~^, then, at least in the case of a sufficiently large initial 
condition, the particle densities no longer approximate a solution of (jl.ip in which the mass 
Io° /m"^ ^fm{x, t)dxdt is conserved throughout time. We refer to [10] and the introduction of [3] 
for a more thorough discussion. 



Our second result shows that the macroscopic density is absolutely continuous with respect to 
Lebesgue measure dn. We will require 

Hypothesis 1.3. There exists a continuous function r : (0, oo) (0, oo) for which T{n)dn = 1, 
with 

/■oo /■ 

/ / (|a;p + I logr(n)| + I log/i„|) /i„ dxdn < oo. 
Jo Jw^ 

and 

(1.8) / / p{n)hn{x)dxdn < oo, 

Jo Jr^ 

where 

n T{m)T{n-m) 
p[n) = / a{m,n — m) --^ dm. 



T{n) 

We also assume that D = sup„ d{'m) < oo. 

Remark. For a simple example for r, consider T{n) = (n+l)~^. If for example a(m, n) < C{m+n), 
then p(n) < Cn and p.8p requires that the total mass to be finite initially. 



Theorem 1.2 Assume that the model is random, and that the law 77V of the initial total particle 
number has a Poisson distribution. Assume also Hypothesis 1.3. Then every limit point V of the 
sequence is concentrated on measures that take the form g{dx,dn,t)dt = fn{x,t)dndxdt, where 
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/ solves (dip. Moreover, there exists a constant C , that may he chosen independently of V, such 
that 



p r /"cxj /• 

(1.9) III '^ifn{x,t))r{x,n) dxdn 



r{dfi) < c, 



for every t, where ip{f) = flog f — f + 1 and r{x,n) = (2tt) '^/^ exp(— |xp/2)r(n). 

Remark. At the expense of discussing some extra technicalities, the proof of Theorem 11.21 might 
include the random model with some other choice of 77V. We only need to assume that for every 
positive A, there exists a constant a(A), such that logE^v exp(A7V) < Na{X). 

Theorem 11.21 is proved by firstly establishing an entropy bound for the distribution of q(t), 
and then using large deviation techniques to deduce that any limit point V of the sequence is 
concentrated on the space of measures g{dx,dn,t)dt = fn{x,t)dxdndt. For this, we simply follow 
the classical work of Guo-Papanicolaou-Varadhan Even though our result is valid for more 
general initial randomness, we prefer to state and prove our results for Poisson-type distributions, 
thereby focussing on the main idea of the method of proof. 

The function r : (0, 00) 1-^ (0, 00) appearing in Hypothesis 1.3 is used to define a reference 
measure with respect to which the corresponding entropy per particle is uniformly finite as e — > 0. 
For simplicity, we take the reference measure which induces a Poisson law of intensity 1 for J\f 
and whose conditional measure i^Ar(-|AA(q) = k) is given by 

k 

(1.10) r(2;j, m-j)(ixi(imj. 

1=1 

The entropy per particle is uniformly finite, because the first part of Hypothesis 1.3 implies that 



supe"'"^ / log F^duN < 00, 
N J 

where F^{q)iyN{dq) denotes the law of q(0). The second part of Hypothesis 1.3 will be used to 
control the time derivative of the entropy. 

We now comment on the possible uniqueness of the solution that the microscopic model ap- 
proximates. We expect to have a unique solution of the system (II. ip for the initial condition h 
as above. However, with the aid of the arguments of [6] and |11] . we know how to establish this 
uniqueness only if we assume that the initial condition satisfies the bound 

POO 

(1.11) / ?T-''||/in||L°° dn < 00, 

Jo 

for sufficiently large b = h{a) (see [6] and [H] for an expression for 6(a)). Using this uniqueness, 
we can assert that in fact the limit V of exists and is concentrated on the single measure 
fj,{dx,dn,dt) = fn{x,t)dxdndt, where / is the unique solution to (jl.ip . As a corollary we have. 
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Corollary 1.1 Assume that Hypotheses 1.1, 1.2 and 1.3 hold and that (jl.lip holds for sufficiently 
large b. Let J : x (0, oo) x [0, oo) —>■ be a bounded continuous function of compact support. 
Then, 

(1.12) limsupETv 

Af^oo 

In (frT2D . f ■.R'^ X [0,oo) X [0,oo) [0,oo) denotes the unique solution to the system (jl.ip with 
the initial data /n(')O) = /i„(-). 

The paper contains an appendix that discusses the scaUngs available in the Smoluchowski 
equations in their continuous form. Examining these scalings produces an heuristic argument for 
the regime of choices of the asymptotic behaviour of the input parameters /3 : (0, oo)^ (0, oo) 
and d : (0, oo) (0, oo) for which a solution (jl.ip will see most of the mass depart from any given 
compact subset of (0, oo) as time becomes high. 

To outline the remainder of the paper: in Section 2, we explain the strategy of the proof, 
giving an alternative overview to that presented in In this section, we also show how the 
microscopic coagulation rate is comparable to the product of densities and may be replaced with 
an expression that is similar to the term Q in p.ip (see Theorem 2.1). The main technical step 
for such a replacement is a regularity property of the coagulation and is stated as Proposition 2.1. 
In Section 2 the proof of Proposition 2.1 is reduced to a collection of bounds that are stated as 
Lemma 2.1. In Section 3, we establish the crucial correlation bound (1.4). In Section 4, the proof of 
Lemma 2.1 is carried out with the aid of the correlation bounds of Section 3. In Section 5, we show 
how the correlation bounds can be used to establish L^-type bounds on the macroscopic densities. 
Sections 6 and 7 are devoted to the proofs of Theorems 1.1 and 1.2 respectively. 

Acknowledgments. We thank James Colliander and Pierre Germain for valuable comments that 
relate to the discussion in the appendix. We also thank an anonymous referee for a number of 
useful suggestions and comments. 

2 An outline of the proof of the main theorem 

Our aim in this section is to outline the proof of the principal result. Theorem 1.1. The overall 
scheme of the proof is the same as that presented in [4j , and the reader may wish to consult Section 
2 of that paper for another overview. 

Our goal is to show that the empirical measures g''{dx, dn,t) converge to f{x,t,dn)dx, where / 
is some measure- valued weak solution of Smoluchowski's equation (II. ip . To this end, we choose a 
smooth test function J : M'^ x (0, oo) x [0, oo) ^ M of compact support and consider the expression 

y(q,t) =e'^-2^ J(xi,mi,i). 
10 



R<* Jo 



J{x, n, t) (/i(dx, dn, dt) — f{x, n, t)dxdndt) 



Evidently, 

y(q(t),t) = J J{x,n,t)g'^{dx,dn,t). 

Note that 

(2.1) y(q(T), T) = y(q(0), 0) + + K^{Y) + A,(y)) {^{t),t)dt + Mr, 
where Mr is a martingale, where the free-motion term KqY equals 

Aoy(q,t) = e'^"^ ^ (i(mi) A^;^ J{xi,mi,t) = J d{n) A^^ J{x,n,t)g''{dx,dn,t). 

and where the collision term AcY is equal to 

(2.2) AcF(q, t) = e'^"^ ^ a{mi,mj)Ve{xi - Xj)J{xi,mi,Xj,mj,t), 
with 14(2;) = e^^y(x/e), and J{xi,mi,Xj,mj,t) given by 

(2.3) ^ — J(xi,mi + nijA) -\ ^ — J(xj,mi + mj,t) — J(xi,mi,t) — J(xj,'mj,t). 

rrii + "ij rrii + rrij 

Our approach is simply to understand which terms dominate in (|2.ip when the initial particle 
number A'^ is high, and, in this way, to see that the equation emerges from considering (j2.1|) in 
the high A'^ limit. Clearly, we expect the last two terms in (|1.6|) . corresponding to free- motion and 
collision, to arise from the terms in ()2.ip in which the operators Aq or Ac act. The time-derivative 
terms in (ll.6p and ()2.ip also naturally correspond. And indeed, the sum of the second and third 
terms on the right-hand side of (2.1) is already expressed in terms of the empirical measure and 
corresponds to the macroscopic expression 



, d{n)AxjJ{x,n,t)f{x,t,dn)dxdt. 

As we will see in Section 6, the term martingale AIt vanishes as e — > 0. The main challenge comes 
from the fourth term on the right-hand side of (2.1), the collision term. How does its counterpart 
in (jl.6p emerge in the limit of high initial particle number? To answer this, we need to understand 
how to express the time-integral of changes to ^(q, t) resulting from all the collisions occurring in 
the microscopic model. To do so, it is natural to introduce the quantity 

/(x,dn;q) = ^ <5-^e(^)<5^,(dn), 

where ^ : M'^ ^ [0, 00) is a smooth function of compact support with ^dx = 1. For 6 > fixed 
and small, in essence counts the number of particles in a small macroscopic region about any 
given point, this region having diamater of order 5. To find the analytic collision term in (|1.6p from 
its microscopic counterpart in ()2.ip . we must approximate the time integral of AcY(q(t),t) by some 
functional of the macroscopically smeared particle count f^, in such a way that the approximation 
becomes good if we take the smearing parameter (5 — > after taking the initial particle number A^ to 
be high. This is achieved by the following important result, in which we write r(q, t) = AcY(ci,t). 
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Theorem 2.1 Assume that the function J(x, m, y, n, t) vanishes when t > T, or m + n < L 
max(m,n) > L. Then 



or 



lim lim sup E ^ 

5^0 AT^oo 



L 



r(q(t),t) -f'^(q(t),t) dt = 0, 



with 



(2.4) 



JK.<i JR<* JO JO 



■oo 



a{m, n)U^ ^{wi—W2)J{wi,m, W2,n, t)f^{wi, dm; q}f^{w2, dn; q)dwidw2 



where we set 



Um,nix) = V{x)[l + u{x;m,n)], U^„^{x)=e '^Um,n{x/e) 



with u{-;m,n) being given in Theorem 1.1. 
Remarks. 

• Note that even thought J is of compact support, the function J given in (2.3) is not in general 
of compact support. In fact, if Xj which appears in (2.2) belongs to the bounded support of J, 
then Xi belongs to a bounded set because of the presence of the term Ve- The same reasoning 
does not work for nii or nij. Of course if J{x,n,t) vanishes if either n > L or n < L~^, then 
J{x, m, y, n, t) vanishes if m + n < L^^, or max(m, n) > L. However, for Theorem 2.1 we assume 
that in fact J vanishes even if one of m or n is larger than L. Because of this, we need to show 
that the contribution of particles with large sizes is small. We leave this issue for Section 6. (See 
Lemma 6.1.) 

• As we mentioned in Section 1, the continuity with respect to m and n and other properties of 
u{-;m,n) will be stated in Lemma 4.2. 

We now explain heuristically why the relation between the cumulative microscopic coagulation 
rate r(q(t),t) and its macroscopically smeared counterpart T^{q{t),t) holds. 

Here is a naive argument that proposes a form for r^(q(t), t). In the microscopic model, particles 
at {wi,m) and {w2,n) are liable to coagulate if their locations differ on the scale of e, \wi — W2\ = 
0{e). If two particles are so located, they coagulate at a Poisson rate of a(m, n)Ve{wi —W2). When 
such a pair does so, it effects a change in y(q, t) of J{wi,m, W2,n). The density for the presence of 
a particle of mass m at location wi should be well approximated by the particle count f^{wi,dm) 
computed on a small macroscopic scale. Multiplying the factors, and integrating over space, we 
seem to show that the expression for r^(q(t),t) should be given by 




where V^{x) = e '^V{x/e). The integrand differs from the correct expression in ()2.4p by the lack 
of a factor of 1 + e'^'^u{{wi — W2)/e;m,n^ . Why is the preceding argument wrong? The reason 
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is the following. The joint density for particle presence (of masses m and n) at wi and 11)2, (with 
\w2 — wi\ = 0(e)) is not well-approximated by the product f^{wi,dm)f^{w2,dn), because some 
positive fraction of particle pairs at displacement of order e do not in fact contribute, since such 
pairs were liable to coagulate in the preceding instants of time, and, had they done so, they would 
no longer exist in the model. The correction factor l+e~'^u[{wi — W2)/e',m,n) measures the fraction 
of pairs of particles, one with diffusion rate d{m), the other, d{n), that survive without coagulating 
to reach a relative displacement wi — 102, and is bounded away from 1 in a neighbourhood of the 
origin of order e. 

We note that in Theorem 2.1 we have reached our main goals, namely we have produced a 
quadratic expression of the densities and a function all which has the macroscopic coagulation 
propensity P for its average. 

The following proposition is the key to proving Theorem 12.11 



Proposition 2.1 Choose T large enough so that J{-,t) = when t >T. We have 



(2.5) 

where 
(2.6) 



lim limsupEjv 



[r(q(t),t)-f,(q(t),t)]dt 



In the statement, z plays the role of a small macroscopic displacement, taken to zero after the limit 
of high initial particle number is taken in the microscopic model. The proposition shows that the 
cumulative influence of coagulations in space and time on Y(q(t),t) is similar to that computed by 
instead considering pairs of particles at the fixed small macroscopic distance z, with a modification 
in the coagulation propensity in the expression ()2.6p being made for the reason just described. 

It is not hard to deduce Theorem 2.1 from Proposition 12.11 We refer to Section 3.5 of ^ for 
a proof of Theorem 2.1 assuming Proposition 12.11 See also [10] for a repetition of this proof and 
more heuristic discussions about the strategy of the proof. 

We will prove Proposition 12.11 in the following way. Define 

^^(q,*) = e^^"^"^^ ^ u^{xi - Xj + z;mi,mj) J{xi,mi,Xj,mj,t), 

where u'^{x;m,n) = e'^~'^u{x/e;m,n). Note that u'^{x) = u'^{x;m,n) solves 
(2.7) (d(m) + d{n))Au^ = a{m, n)(Ku^ + V^), 

with 

V,{x) = e-V(x/e), V'{x) = e-'^Vix/e). 
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The process { {X^ - Xq) (q(t), i) : t > 0} satisfies 

(2.8) (X, -Xo)(q(T),r) = (X,-Xo)(q(0),0) (- + Aq) (X, - Xq) (q(t), 



+ / A,{X,-Xo){q{t),t)dt + M(r), 
Jo 



with {M(t) : t > 0} being a martingale. We will see that the form ()2.5p emerges from the dominant 
terms in ()2.8p . those that remain after the limit of high initial particle number — 00 is taken. 
To see this, we label the various terms which appear on the right-hand side of ()2.8p . Firstly, those 
terms arising from the action of the diffusion operator: 



d_ 

dt 



Ao {X, - Xo) = Hii + H12 + + Hu + + + 



with 
Hnici,t) 



^2{d 2) a{mi,mj) [xi — Xj + z) — (^Xi — xj) J{xi,mi,Xj,mj,t), 
_^2{d-2) a{mi,mj)V):[xi — Xj)u'^{xi — Xj;mi,mj^J{xi,mi, Xj,mj,t), 



a{mi,mj)V^{xi — Xj + — Xj + z; mi,mj)j{xi, nii, Xj,mj,t), 



-f^i4(q)*) = e^^'^ d{mi) u'^{xi — Xj + z;mi^mj) — u'^{xi — Xj;mi,mj) Jt{xi,'mi,Xj,'mj,t), 
along with 

H2{^.,t) = 2e^^'^~'^^ d{mi) u%{xi — Xj + z]mi,mj) — u%{xi — Xj^nii^rrij) ■ Jx{xi,mi,Xj,mj,t) 

H3{(i,t) = — 2e^*-'^^^^ d{mj) ul.{xi — xj + z;mi,mj) — ul.{xi — Xj;mi,mj) ■ Jy{xi,mi,Xj,mj, 



and 



-f^4(q) = e^^'' — Xj + z; mijUij) — u'^{xi — Xj] mi,mj) 



«jG/q 



\^d{mi)AxJ{xi, mi,Xj,mj,t) + d{mj)/S.yJ{xi, mi,Xj,mj,t) 



where Jx denotes the gradient of J with respect to its first spatial argument, Jy the gradient of J 
with respect to its second spatial argument, and • the scalar product. As for those terms arising 
from the action of the collision operator. 



Ae(X, - Xo){ci, t)=Gl (q, t) +Gl{q,t)- Gl (q, t) - G^ (q, t) , 
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where ^^(q, t) is set equal to 

ruk 



rrik + mi 



+ 



mi 



irik + mi 



u^{xk - Xi + z;mk + mi, mi)J{xk,mk + mi, Xi,mi,t) 

+ - Xk + z;mi,mk + me)J{xi,mi,Xk,mk + mi,t) 
u^xi - Xi + z;mk+ mi, mi)J{xi, nik + mi, Xi,mi,t) 

+ u''{xi - Xi + z;mi,mk + mi)J{xi,mi,xi,mk + mi,i) 
u'^ixk - Xi + z; mk,mi)J{xk,mk,Xi,mi, t) 

+ - Xk + z;mi,mk)J{xi,mi,Xk,mk,t) 
u^ixi - Xi + z; rrii, mi)J{xi, nii, Xi,mi,t) 

+ u''{xi - Xi + z;mi,mi)J{xi,mi,Xi,mi,t) 



and where 

G1{(\, t) = -e^('^~^) ^ a{mk,mi)Ve{xk - xe)u^{xk - xi + z; mk,mi)J {xk,mk,xi, mi, t). 

The terms in G\ arise from the changes in the functional Xz when a collision occurs due to the 
influence of the appearance and disppearance of particles on other particles that are not directly 
involved. Those in are due to the absence after collision of the summand in indexed by the 
colliding particles. 

As we take a high N limit in (|2.8|) . note that the quantity 
rT 

/ r(q(t), = e^^'^"^-* a{mi,mj)V'^{xi — Xj)J{xi,mi,Xj,mj,t) 

appears, with a negative sign, in the term Hn. The term H12 also remains of unit order in the high 
limit, and would disrupt our aim of approximating T{ci{t),t)dt by z-displayed expressions. 
However, our definition of n*^ (see (j'2.7|) ) ensures that 

H12 — Gq = 0, 

so that this unwanted term disappears. The definition of u'' was made in order to achieve this. The 
other term of unit order remaining in the high N limit is the z-displaced Hi^. Rearranging ()2.8p . 
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we obtain 



/ Hn{cL{t),t)dt+ [ HMt),t)dt < |X, -Xo|(q(r),T) + !X, -Xo!(q(0),0) 
Jo Jo 

(2.9) +/ {\Hu\ + \H2\ + \H3\ + \H4){ci{t),t)dt 

Jo 

+ [ \Gl-Gl\{qit),t)dt+ [ \GlMt),t)dt 
Jo Jo 

+ \M{T)\. 

We have succeeded in writing T^ — F in the form Hn+Hi^, so that, for Proposition l2.H it remains to 
prove that the right-hand-side of ()2.9p is small enough. Firstly, recall that, by our assumption, the 
function J is of compact support. We now choose T sufficiently large so that J{x, m, y, n, T) = 0. 
As a result, the first term on the right-hand side vanishes. The other bounds we require are now 
stated. 

Lemma 2.1 There exists a constant C2 = C2{J,T) such that, 

(2.10) (IF2I + \H3\) {il{t),t)dt < C2\z\^, 

Jo 

(2.11) / EN{\H^\ + \Hu\){ci{t),t)dt<C2\z\^\ 

Jo 

(2.12) / EN\Gl-GlMt),t)dt<C2\z\^, 

Jo 

(2.13) / E\Gl\{q{t),t)dt<C2(^^ 



JO ^ I- 

(2.14) EMX, - Xo\{qm < C2\z\, 

(2.15) Eiv[M(T)2] < C2e'^-^. 

These bounds are furnished by the correlation inequality Theorem 13.11 that is the main innovation 
of this paper, to whose proof we now turn. 



3 Correlation Bounds 

This section is devoted to the proof of the correlation bound which appeared as (1.4) when d{-) 
is non-increasing and takes the form (3.1) in general. Recall the function (p which appeared in 
Hypothesis 1.1. The main result of this section is Theorem 13.11 
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Theorem 3.1 For every non-negative bounded continuous function K : 
(3.1) E;v / K{xi,{t),...,x,^{t))\{-fk{mat))dt 

° n,...,«fcG/q(t) r=l 

k 

<En (^"'^ ^°^^) (^n (0) , . . . , X,, (0)) n Ik (m,, (0)) , 

jl,...,ifcG/q(o) r-=l 

where all summations are over distinct indices ii, . . . ,i)^, the function 7fc(?TT.) = md{m)'^/'^(f){m)~~^ , 

and the operator A is defined by 

(3.2) 

A-^'-'«'=i^(yi,...,yfc) = co(A;d) | J] d(n,)-'^/2 (^^^^-^ + • • • + ^^^-i^) ' i^(zi, . . . , z,)dz„ 
where Co{kd) = {kd — 2)~^uj~f^^ , with uJkd denoting the surface area of the unit sphere in 



Let us make a comment about the form of ()3.ip before embarking on its proof. Observe that 
if there were no coagulation, then it would have been straightforward to bound the left-hand side 
of ()3.ip with the aid of the diffusion semigroup even if we allow a function K that depends on 
the masses of particles. Indeed, if ' denotes the diffusion semigroup associated with 

particles (xj^,mj^), . . . , {xi^,^mi^), then j'^ 5™'^' is exactly the operator A™'i''"'™''fe . What 
(3.1) asserts is that a similar bound is valid in spite of coagulation provided that we allow only a 
very special dependence on the masses of particles. 

Proof of Theorem 3.1. Let us define 

k 

G(q)= Y (A"^--'"^^^i^)(x.„...,x,jn^'«("^v)- 
n,...,ifce/q r=i 

Recall that the process q(t) is a Markov process with generator L = Aq + Ac where Aq = 
Yjiei^d{'mi)l^x^. By Semigroup Theory, 

(3.3) EivG(q(t)) =E7vG'(q(0)) / LG(q(t))dt. 
We have 

k 

(3.4) ^oG(q) = - Y K{xi^,...,Xi^)J{-ik{mi^). 

ii,...,«fce/q r=l 

This and the assumption K >Q would imply (j3.ip provided that we can show 

(3.5) AcG < 0. 

To prove (3.5), let us study the effect of a coagulation between the i-th and j-th particle on G. 
We need to study three cases separately: 
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• hj ^ {h,---,ik}, 

• i,j G {ii,...,ik}, 

• only one of i or j belongs to {ii, . . . ,ik}- 



If the first case occurs, then (i, j)-coagulation does not affect the term indexed by {ii,. . . 
in G(q). 

If the second case occurs, then we need to remove those terms in the summation for which 
{hj} C {ii, . . . , ifc}. This contributes negatively to AcG(q), because K >0. This total contribution 
for this case is given by 

- X] Ve{xi - Xj)a{mi,mj) 

k 

■ e{il,...,ik}){A"''^'■■■'"''^K){xi,,...,Xi^)]JJk{rni^)■ 
^l,■■■,^k r=l 

If the third case occurs, then only one of i,j belongs to {ii, . . . , i^}. For example, either i = ii, 
and j ^ {ii, . . . , i^}, or j = h, and i ^ {ii, . . . , i^}- In this case, the contribution is still non-positive 
because after the coagulation the expression 

k 

k 

i2,..;ik r=2 

is replaced with the expression Y2 which is given by 

k 



■> 1-2, ■■■,tk r=2 

k 

+ i2,...,ik r=2 

For (3.5), it suffices to show that Y2 < Yi. For this, it suffices to show that for every positive 
m, n, A and B, 



(3.6) ^(m + n)^-^ 



-| kd 

a[m + n) 



, . . kd 1 r . 1 kd 

< (i){m)^-^[A + B]^-^. 



We are done because the assertion (3.6) for fixed m, n and all positive A and B is equivalent to the 
inequalities 

(t){m)d{m) > (l)(m + n)d(m + n), 
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and 

<f>{m) > (j){m + n), 

both being satisfied, and these are true for all choices of m and n by Hypothesis 1.1. □ 
Corollary 3.1 For every non-negative bounded continuous function K, 



(3.7) s'^me^t ^ Yl ^KW' ■ ■ -^^iM n 

^ ii,.--,ifee/q(t) r=l 

< Co (fed) / K{xi, ...,Xk)Y[{hk* Afe) {Xr)dXr, 
J r=l 

where hk = n^{n)~2 ^d{n)^~khn dn and \k{w) = \w\k'''^. 



Proof. Prom the elementary inequality ai . . . < (af + • • • + o^)'^''^, we deduce that the kernel 
^ni,...,nfc q£ ^Yie operator A^i' - '"*^ is bounded above by 

k 

A"!'-'"^-^!, ■■■,zk)< co{kd) H \zr\-^-'^d{nr)-^ . 

r=l 

This and (3.1) imply (3.7). □ 
We end this section with two lemmas concerning the first condition in Hypothesis 1.1. 

Lemma 3.1 Suppose the function d{-) has a finite negative variation in an interval [a,b] C (0, oo). 
Then there exists a positive continuous function (p such that cp and (pd are non-increasing in the 
interval [a,b]. 

Proof. Step 1. Firstly, we assume that there exist points ao = b > ai > ■ ■ ■ > a^-i > a( = a 
such that d(-) is monotone on each interval [aj,ai_i], i = For the sake of definiteness, 

let us assume that d{-) is non-decreasing (non-increasing) in [aj,aj_i], if i is odd (even). In this 
case, we can construct a continuous as follows: Define Aq = A and A^. = ^ITiLi dfi^t^'i) 
k > 1. For X G [o2A;+i, a2A;] and A; > 0, we set 0(x) = For x G [02A;) a2A;-i] and A; > 1, we set 

Step 2. Let d be a continuous positive function. Approximate d in L°° by a sequence of continuous 
piecewise monotone functions {dn}. To simplify the presentation, we assume that each dn is as in 
Step 1. That is, dn increases near the end point b. Let us write cpn for the corresponding cp, and let 
Cn denote the number of intervals in the partition (so that ac„ = a). It remains to show that the 
sequence {cpn} has a convergent subsequence. Since each (pn is non-increasing, we may appeal to 
the Helley Selection Theorem. For this we need to make sure that the sequence {(pn} is bounded. 
Note that supj.g[„;,] (pn{x) = (pn{o-) = (pn{o'c„)- Set Dn = A cn-i if Cn is odd and Dn = ^^-1 if 
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Cn is even. We readily see that (/>„(c„) < ( inf^g[a ^ (i(x)) L>„, whatever the parity of Cn- The 
infimum being positive, we require that sup„gj^Z)„ < oo. For any A; G N for which Ak is defined, 
we may take the logarithm of to produce a sum and observe that d(-) is non- increasing on the 
intervals [a2i,a2i-i]- Hence, log^^ measures the negative variation of the function logd on the 
interval [a2k, b\ . Since d is uniformly positive, sup„ Dn < oo is implied by the function d having a 
finite negative variation. □ 

Lemma 3.2 Suppose the function logd{-) has a finite negative variation in an interval [no, oo) with 
no > 0. Then there exists a function positive continuous (j) such that (j) and (j)d are non-increasing 
in the interval [no,oo). 

Proof. The proof is very similar to the proof of Lemma 3.1. First we assume that d is piecewise 
monotone. This time we set (pino) = A and define (j) continuously so that (f) is constant when 
d decreases and (/> is a constant multiple of d~^ when d increases. Since is non-increasing, we 
may end with a function which crosses and becomes negative. This can be fixed by adjusting 
A = (j){no), only if (p is bounded below. As in the proof of Lemma 3.1, we can readily see that (p is 
bounded below if the total negative variation of log d is finite. □ 

Note that in the statement of Lemma 3.2 we can not drop log because on the infinite interval 
[nj,oo) the function d(-) could take arbitrarily small values. 



4 Proof of Lemma 2.1 



The strategy of the proof of Lemma 2.1 is the same as the one used to prove the analogous 
inequalities in [3]. The only difference is that we only need to use our correlation bound Corollary 
3.1 to get the bounds (2.10-15). For (2.10) and (2.11) we need to apply Corollary 3.1 for k = 2. 
Corollary 3.1 in the case k = 3 will be used for (2.12). As for (2.15) all cases k = 2,3,4 will be 
employed. We omit the proof of the inequalities (2.13) and (2.14) because they can be established 
by a verbatim argument as in In fact the proof (2.14) is straightforward because we are dealing 
with a calculation involving the initial configuration. For this, however, a suitable bound on the 
function would be needed that will be stated as a part of Lemma 4.2 below. The same bound 
and Lemma 4.1 below will imply (2.13). 

The main ingredients for the proof of inequalities (2.10) and (2.11) are Corollary 3.1 (with 
k = 2), certain bounds on and u% (which will appear in Lemma 4.2), and Lemma 4.1 below. 
The straightforward proof of Lemma 4. 1 is also omited and can be proved in exactly the same way 
we proved Lemma 3.1 of 0]. 

Lemma 4.1 For any T € [0, oo), 




As for the remaining inequalities, we only establish (2.12) and (2.15) because these are the most 
technically involved cases and the same idea of proof applies to (2.10) and (2.11). 

We now state our lemma about the functions u and n^. Recall that u'^{x; n, m) = e^~'^u{x j e\ n, m) 
where u satisfies 

Au(x; n, m) = a'{n, m)V{x) [l + u{x; n, m)] , 
with u{x;n,m) ^ as |x| oo, and 

, a{n,m) 
a [n,m) := — — -. 

For our purposes, let us write for the unique solution of 

Aw^ix) =aV{x)[l + w''{x)], 

with w^{x) ^ as |x| — > oo. Of course, if we choose a = a'{n,m), then we obtain u{x;n,m). We 
choose the constant Co so that V{x) = whenever |x| > Cq. 

Lemma 4.2 There exists a constant C3 for which the following bounds hold. 

• — 1 < w°'{x) < and for x € M"^, 

k"(x)| < C73amin{|x[2-'^,1}, 
K(^)l < C3amin{|x|^"'^,l}. 

• for x € satisfying \x\ > max|2[z| + Coe, 2Coe}, 

(4.1) |?x'^(a; + z; n, m) — u'^{x; n,m)\ < C^a {n., m)\z\ |x|^~'^ 
and 

(4.2) |u|(x + n, m) — n|(x; n, m)| < C^a' {n,m)\z\\x\^'^ . 

• the function w"" is differentiable with respect to a and a~^w°' < < 0. 

Proof. The proof of the first and second parts can be found in Section 3.2 of [1| and we do not 
repeat it here. As for the third part, recall that the function w"" is uniquely determined by the 
equation 

(4.3) w^{x) = -coa [ \x- 2/p-''F(y)(l + w''{y))dy, 

where cq = co{d) = {d — 2)^^uj^^, with tOd denoting the surface area of the unit sphere S'^^^. We 
wish to show the regularity of the function w"" with respect to the variable a. In fact the existence 
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of the unique solution to (4.3) was established in [Ij using the Fredholm Alternative Theorem. To 
explain this, let us pick a bounded continuous function R such that R > 0, with 



/ R{x)dx = oo, / R{x)\x\'^ '^'^dx<oo. 



Define 



7^=<n:M'^^M: u is measurable and / {x)R{x)dx < cx3 > . 
Observe that 7^ is a Hilbert space with respect to the inner product 

{u,v) = / u{x)v{x)R{x)dx. 

Note that if w'^ solves (4.3), then, defining T -.Ti ^ Tihy 

= co [\x- y\^-''V{yMy)dy, 



we have that 

(4.4) (id + a.F)(u;°) = -ar 

where 

r(x) = CO / \x- y\^-''V{y)dy, 



and id means the identity transformation. We wish to show the differetiability of w"' with respect 

da ' 



to a > 0. This is clear heuristically because we have a candidate for v"" := if we differentiate 



both sides of (4.4), then v"" solves 

(4.5) {id + aT){v^) = -F - Tw"" = a'^w''. 

This provides us with a candidate for , because the operator id + aJ-" has a bounded inverse 
(see Section 6 of [4J). The rigorous proof of the differentiability of w"" goes as follows. First define 
^a,/i _ (^^a+h _ ^ayj^ observe that v""'^ satisfies 

(4.6) {id + aT){v'''^) = -T - Tw^^^. 

We would like to show that v""'^ has a limit in 7i, as h ^ 0. One can readily show that the 
right-hand side of (4.6) is bounded in 7i because |i(;'^(x)| < C2amin{|xp~'^, 1} by the first part of 
the lemma. Hence v""'^ stays bounded as /i — > 0. If is any weak limit, then v"" must satisfy 
(4.5). Since (4.5) has a unique solution, the weak limit of v""'^ exists. In [3], it is shown that is a 
compact operator. From this and (4.6), we can readily deduce that the strong limit of v""'^ exists. 
As a consequence, w"" is weakly differentiable in a and its derivative satisfies (4.5). Using Sobolev's 
inequalities and the fact that V is Holder continuous, we can deduce by standard arguments that 
indeed is and satisfies 

(4.7) Av^ = av^V + (1 + w'')V. 
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This means that w"'{x) is continuously differentiable with respect to (x,a). 

We now want to use (j4.7p or equivalently (|4.5p to conclude that a~^'w"' < v"" < 0. In fact, by 
■Sp . we have that = —aTv^ — a~^w"', which implies that 

2-dj„. I „-l|„.,a/ 



< c^co y |x-yr"'^(iy + a"V(2;)l, 

where c'^ is an upper bound for with x in the support of the function V. From this, it is not 

hard to deduce that there exists a constant such that 

(4.8) |t;"(x)| < max{\x\'^-'^ , 1}. 

In a similar fashion, we can show that there exists a constant such that 

(4.9) \Vv''{x)\ < cl maxjlxl^""', 1}. 

We now demonstrate that v"" < 0. Take a smooth function 995 : M ^ [0, 00) such that ip'g,fs>0 
and 

fo r<0. 



r r > 6. 



We then have 

(4.10) - / ip's{v'')\Vv''\'^dx = [ ipsiv'')Av''dx = [ Vil + w" + av'')ip5ivldx, 

Jr^ jRd J^d 

the second equality by (j4.7p . Integration by parts was performed in the first inequality: we write 
the analogue of (|4.10p which is integrated over a bounded set {x : \x\ < R} . We may obtain (j4.10p 
by sending R ^ 00 but for this we need to make sure that the boundary contribution coming from 
the set {x : \x\ = R} goes away as ^ 00. This is readily achieved with the aid of (j4.9p . Since 
1 + tL)" > by the first part of the lemma, and v°'(ps{v°') > 0, we deduce that the right-hand side 
of (I4.10p is non-negative. Since the left-hand side is non-positive, we deduce that 



/ 



ip's{v'')\Vv''\'^dx= I Vil + w" + av'')(ps{v'')dx = 0. 
We now send (5 — > to deduce 



0=/ iVv^l^liv" >0)dx = I Vil + w" + av'')v''l{v'' >0)dx. 

As a result, on the set A = {x : v"" > 0} we have Vf " = 0. Hence v"" is constant on each component 
B of A. But this constant can only be because on the boundary of A we have v"" = 0. This is 
impossible unless A is empty. Hence, < everywhere. 

It remains to prove that > a^^w'^. For this observe that if 7° = a~^'w'^ — v"", then 



= aV-i" + V{-w''). 
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We can now repeat the proof of < to deduce that 7" < because —w°' > 0. This completes 
the proof of the third part of the lemma. □ 

Proof of (2.12). Note that 



I ^N\Gl-Gl\{ci{t),t)dt<y^Di, 
•^0 i=i 



where the first four of the Di are given by 



rrik + me. 



Di = En dt a{mk,'me)Ve{xk - ■ 

{u^ixk -Xi + z;mk + mi,mi) - - xfjirik + me,mi)\ \j{xk,mk + rrii, Xi,mi,t)\, 

[u'^ixe -Xi + z;mk + me,mi) - u^ixi - Xi\mk + m£,mi)| \ j{xi,mk + mi,Xi,mi,t)\, 

i€/q 

rT 

D3 = En dt y2 a{mk,me)V,{xk-xe)e^^^-^^ 



{u^ixk - Xi + z;mk,mi) - u^{xk - Xi;mk,mi)\ \ j{xk,mk, Xi,mi,t)\, 

ieiq 



and 

rT 



D4 = En [ dt y2 a{mk,me)Ve{xk - xi)e^^'^-^^ 



fc,^e/q 



y^ \u^ixe - Xi + z;mt,mi) - u^{xt - Xi;me,mi)\ \j{xe,me,Xi,mi)\. 



The other four terms each take the form of one of the above terms, the particles indices that appear 
in the arguments of the functions and J being switched, along with the mass pair labels for these 
functions. 

The estimates involved for each of the eight cases are in essence identical. We will examine 
the case of D^. We write Ds = + D"^, decomposing the inner i-indexed sum according to the 
respective index sets 

{i G /q, z 7^ k,i, \xk — Xi\ > p} andji E Iq,i k,£, \xk — Xi\ < p} 

Here, p is a positive parameter that satisfies the bound p > max {2\z\ + Cqc, 2Coe}. By the second 
part of Lemma 4.2, we have that 

rT 

' xe) 



< ° — Eat dt a{mk,mi)V^{xk 

P Jo u 6^T 



k,eeiq 
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where we have also used the fact that the test function J is of compact support, and the fact that 
the total number of particles living at any given time is bounded above by Ze^~'^. Prom the bound 
on the collision that is provided by Lemma 4.1, follows 

^1 ci\z\ 



1 ■ 



To bound the term D^, note that by Lemma 4.2, the term is bounded above by 

• ^ Ti-ilxi - Xk\ < p} ^[u^ixk - Xi + z;mk,mi)\ + [u^ixk - Xi;mk,mi)\ \j{xi,mi,Xk,mk,t)\dt 
<ciEiv / e3('^-2) (ximk,mi)V'{xk-Xi) 

• Y ~ < P, max {mfc,mi, \xk\, \xi\} < L, nik + rrii > L~^} 



a'{mk,mi) 



I'^k Xi -\- Z\ -\- \Xj^ Xi\ 



dt. 



where = e^^'^Ve and L is chosen so that J(a;, m, y, n) = if any of the conditions 

m + n>L~^, inax{m,n) < L, max(|x|, |y|) < L, 

does not hold. We note that if m + k + nii > L~^, then a'{mk, rrii) < C2a{mk, rrii), for a constant 
C2 that depends on L. On the other hand, the conditions 

rrik < L, rrii < L, rrik or rui > 2^~^ i 
imply that for a constant C3 = C3(L), 

a(mk,mi)a{m,k,mi) < C3^3{mih3{mi)'~f3{mk), 

where we have used second part of Hypothesis 1.1. We are now in a position to apply Corollary 
3.1. For this we choose k = 3 and 



K{xi,X2,X3) = V^{xi - X2)1{\X2 - X3\ < p, \X2\, \X3\ < L} 

As a result, < D{z) + D{0) where D{z) is given by 



I , |2-d , I \2-d 

\X2—X3 + Z\ +\X2—X3\ 



f 2-d ^ - 

C4 {Xi - X2)1{\X2 - X3\ < p, \X2\,\X3\ < L}\X2 - X3 + z\ JJ (/i.3 * A3) (Xr)^: 

1 

< C5 y V^{xi - X2)1{\X2 - X3I < p, \X2\, \X3\ < L}\ 



X2 — X3 + z\'^ '^dxidx2dx3 



< C6 / la + zp '^da < C7{p+ \z\f, 

J\a\<p 



25 



where, for the first inequahty, we used Hypothesis 1.2(ii). Combining these estimates yields 

D^ = D^ + D^ <ci^ + cr{p+\z\)\ 

1 2 

Making the choice p = \z\'>-+'^ leads to the inequality < csl-zl'^+i. Since each of the cases of 
jz^j : i G {1, . . . , 8}| may be treated by a nearly verbatim proof, we are done. □ 

Proof of (2.15). Setting L = Aq + Ac, the process 

M,(r) = x,(q(r),r)-x,(q(o),o)- (^- + I.)x,{q{t),t)dt 

is a martingale which satisfies 

En [M,{Tf] = En C (LX^ - 2XJuX,) {(i{t),t)dt = V E,v A(q(t), t)dt, 
Jo .^^ Jo 

where 



Aiiq,t) = 2e^^'^ ^) ^ d{mi) V xi ^ u''{xi - Xj + z]mi,mj)J{xi,mi,Xj,mj,t) 



and 



^2(q,i) = 2e^('^ ^ d(mj) V^, ^ - + z;mi,mj)J{xi,mi,Xj,mj,t) 
while ^3(q, t) is given by 



(4.11) 



4(d-2) 



a{mi,mj)e "^Veixi-Xj 



rrii 



rrii + ruj 



u^{xi -Xk + z;mi + mj,mk)J{xi,mi + mj,Xk,mk,t) 



rrii 



+ 



rrii + nT-j 
rrii 



rrii + ITT'- 



+- 



m 



■3 



u^^Xk -Xi + z;mk,mi + mj)J{xk,mk,Xi,mi + mj,t) 

-u''{xj - Xk + z;mi + mj,mk)J{xj,mi + mj,Xk,mk,t) 
i 

-u^{xk — Xj + z; rrik, rrii + fnj)J{x}~, m^, Xj,mi + mj,t) 



rrii + 

-u^ixi - Xk + z;mi,mk)Jixi,mi,Xk,mk,t) 
-u'^ixk - Xi + z;mk, rrii) J{xk,mk,Xi, rrii, t) 
-u^ixj - Xk + z;mj,mk)J{xj,mj,Xk,mk,t) 
-u^{xk-x, + z;mk,m,)J{xk,mk,x„m„t)^ 

—u^{xi — Xj + z; rrii, mj)J{xi, rrii, Xj, rrij, t)| 
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We now bound the three terms. Of the first two, wc treat only Ai, the other being bounded by an 
identical argument. By multiplying out the brackets appearing in the definition of Ai, and using 
supj„g(o,oo) d{m) < oo, (which is assumed by Hypothesis 1.1), we obtain that Ai < Au + ^12 with 

All = coe^^*^"^^ ^ \ul{xi-Xj + z;mi,mj)\\ul{xi-Xk + z;mi,mk)\ 

■\J{xi,mi,Xj,mj,t)\\J{xi,mi,Xk,'rnk,t)\ 
A12 = coe^^*^"^^ ^ {xi- Xj + z;mi,mj)\\u^ {xi - Xk + z;mi,mk)\ 

•\Jx(^Xi^ TJli , Xj , 771 j , t) 1 1 e/j; i^Xj,^ "^ii Xk J '^kj ^) | • 

Let us assume that z = because this will not affect our arguments. We bound the term An with 
the aid of Corollary 3.1 and Lemma 4.2. The term A12 can be treated likewise. To bound ^n, 
first observe even though i and j are distinct, k and j can coincide. Because of this, let us write 
All = Am + ^112 where Am represents the case of distinct i,j and k. We only show how to 
bound ^111 where the correlation bound in the case of = 3 is used. The term ^112 can be treated 
in the similar fashion with the aid of Corollary 3.1 when k = 2. Since J{x, m, y,n) 7^ implies that 
m, n, \x\, \y\ < L and m + n > L~^. Using second part of Hypothesis 1.1, we can find a constant 
ci = ci(L) such that 

a{7ni,77lj)a{r7li,77lk) < C2^3{mh3iinj)^3{77lk), 

whenever 

rui, mj,mk < L, rrii + mj,mi + rrik > L~^. 
As a result, we may apply Corollary 3.1 with k = 3 and 

K{xi,X2,X3) = e'^~'^\xi - a;2|^"'^|xi - a;3|^"'^ll(|xi|, \x2\, \xs\ < L), 

to deduce 

Am < C2e'^ M |xi - X2\^ '^\xi - X2,\^ '^H\xil \x2\, \x2,\ <L)W^ (/is * A3) {Xr)dXr. 

J r=l 

Note that K is an unbounded function and Corollary 3.1 can not be applied directly. However we 
can approximate K with a sequence of bounded functions and pass to the limit. Prom this and 
Hypothesis 1.2, we deduce 

This and an analogous argument that treats the terms A112, A12 and A2 lead to the conclusion 
that 

(4.12) Ai+A2<cie'^-'^. 

We must treat the third term, ^3. An application of the inequality 

(ai + . . . + Qnf < n{al + . . . + a^) 
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to ^3, given in ()4.1ip . implies that 

9 



3i) 

i,j&Ici n=l k&Iq i=l 



(4.13) ^3(q,t)<96^('^-') J2 c.{m,m,)V,{x,-x,)[^{^Yn)'+Yi] =:^A 

where Yi is given by 



u^{xi - Xk + z;mi + mj,mk)J{xi,mi + mj,Xk,mk,t), 



nii+m 



and where {Yi : i € {2, . . . , 8}} denote the other seven expressions in (j4.1ip that appear in a sum 
over /c G /g, while Ig denotes the last term in (j4.1ip that does not appear in this sum. There are 
nine cases to consider. The first eight are practically identical, and we treat only the fifth. Let us 
again assume that z = because this will not affect our arguments. Note that 



^35 = £^('^-2) J] a(m„m,)V;(x,-x,)(^y5)' 
= e^^'^-^) a{mi,mj)V'{xi-Xj) 

(xi - Xk;mi,mk)u'^ (xj - xi;mi,mi) J{xi,mi,Xk,mk,t)J{xi,mi,xi,mi,t) 

In the sum with indices involving k,l G /q, we permit the possibility that these two may be equal, 
though they must be distinct from each of i and j (which of course must themselves be distinct 
by the overall convention). Let us write ^35 = + A352, where A351 corresponds to the case 
when all the indices i,j,k and I are distinct and ^352 corresponds to the remaining cases. Again, 
our assumption on a as in Hypothesis 1.2 would allow us to treat the term ^4351 with the aid of 
Corollary 3.1. This time k = A and our bound on u given in the first part of Lemma 4.2 suggests 
the following choice for K: 

K{xi, ... ,X4) = e°'"^y^(xi - X2)\xi - X3p~°'|xi - a;4p~'^l(|xi|, \x2\, {x^l, \x^\ < L). 

Note that K is an unbounded function and Corollary 3.1 can not be applied directly. However we 
can approximate K with a sequence of bounded functions and pass to the limit. From Corollary 3.1 
and Hypothesis 1.1 on the initial data we deduce that the expression A^^idt is bounded above 
by 

Cse"^"^ j V^{xi — X2)\xi — X3p~'^|xi — X4p^'^l(|xi|, |X2|, |X3|, |X4| < L)dxi . . . dx4 = CqE 



d-2 



A similar reasoning applies to ^352, except that Corollary 3.1 in the case of /c = 3 would be 
employed. Hence, 



(4.14) J]^3i<C7£'^-' 

i=l 
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We now treat the ninth term, as they are classified in (I4.13p . It takes the form 

g4d-8 a{mi,mj)Ve{xi — Xj)u^{xi — xj + z;mi,mj)'^J{xi,mi,Xj,mj,t)'^. 

This is bounded above by 

cge^'^'^ ^ a{mi,mj)Veixi - Xj), 

because u"" < cge^"'^ by the first part of Lemma 4.2. The expected value of the integral on the 
interval of time [0, T] of this last expression is bounded above by 

fT 

cje^'^~'^EN / y2 a{mi,mj)Ve{xi - Xj)dt < cioe'^~'^. 

where we used Lemma 4.1 for the last inequality. This, (4.12), (4.13) and (4.14) complete the proof 
of (2.15). □ 



5 Bounds on the Macroscopic Densities 

In this section we show how Corollary 3.1 can be used to obtain certain bounds on the macroscopic 
densities. These bounds will be used for the derivation of the macroscopic equation. Recall that 

g^{dx,dn,t) = e'^~^^6(^^^(^t),mi{t)){dx,dn), 

i 

and that the law of 

q I— > g'^{dx, dn, t) 

induces a probability measure on the space X. The main result is Theorem 5.1. 

Theorem 5.1 Let V be a limit point ofV^. The following statements are true: 
• 1. For every positive Li, and k G {2,3,4}, 

dxdtdV < oo, 



5 JXJO J\x\<L 



^1 



oo 



(5.1) sup I I I I U\x-y)jkin)gidy,dn,tf^ 



where ^^{x) = 5 '^^(^j^, with ^ a nonnegative smooth function of compact support satisfying 

2. We have g{dx,dn,t) = f{x,t,dn)dx for almost all g with respect to the probability measure 
V. 
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3. For every continuous R of compact support and positive L, 

(5.2) liinj J J J J R{x,m,n,t)f{x,t,dm)f{x,t,dn)dxdt 



T rL rL 



10 JL-^ JL- 

where 



R{x, m, n, t)f{x, t, dm)f{x, t, dn)dxdt 



dV = 0, 



(5.3) fix, t, dn) = j y)gidy, dn, t). 

Proof. Fix X and choose 

k 

K{yi, . . . , yfc) = JJ - Vr), 

r=l 

in Corollary 3.1. The right-hand side of (3.7) equals 

J H ^^{^ - Xr)hk * \k{Xr)dXr, 



r=l 



which, by the second part of Hypothesis 1.2, is bounded by a constant ci(Li) when A; = 2, 3, 4, and 
|x| < Li. As a result, 

h 

(5.4) Ejv / / €^=('^-2) V r\^'ix-Xi^it)hkimiAt))dxdt<c,{Li) 

for a constant ci(L) which is independent of S and e. Here we are assuming that the indices ii, . . . ,ik 
are distinct. Note that if wc allow non-distinct indices in the summation, then the difference would 
go to as e ^ because the summation is multiplied by £^i^~'^) while the number of additional 
terms is of order 0{e^'^~^^^'^~'^'>). As a consequence, we can use (5.4) to deduce (5.1). 

Recall that the function 7^ is a positive continuous function. From this and (5.1), one can 
readily deduce part 2. 

It remanis to establish part 3. First observe that by (5.1) and the posivity of 74, 



(5.5) sup 



T 



L 



f{x,t,dn) 



5 J Jo J\x\<Li 

Because of this, it suffices to prove that 



4 



dxdtV{dg) < 02(^1, L). 



lim / / / / Rpix, m, n, t)f{x, t, dm)f{x, t, dn)dxdt 

S^oJq J J^-iJ^-i 

= / / Rp{x,'m,n,t)f{x,t,dm)f{x,t,dn)dx. 

Jo Jl-^ Jl-^ J 
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for each p, provided that liuip^oc Rp{x,rn,n,t) = R{x,m,n,t), uniformly for 111,71 ^ [L~^,L], 
\x\ < Li and t < T. By approximation, we may assume that R is of the form R{x,m,n,t) = 
"^1=1 Ji{x,t)J2{m)J^{n). Hence it suffices to estabhsh (\5.2h for R of the form R{x,m,n,t) = 
Ji{x,t)J2{rn)J'i{n). But now the left-hand side of (15. 2|) equals 



lim 

5^0 



Ji(x, t)dxdt. 







L 

J2{fn)f^{x, t, dm) 



L 

Mn)f\x,t,dn) 



We note that ^ ^ 

^J2im)fix,t,dm) = {^j^ ^J2im)fi-,t,dm)^ *xC\x). 

converges almost everywhere to 

J2(m)/(x, t, dm). 

The same comment applies to J^{n)f^{x,t)dn. From this and (5.5) we deduce (5.2). □ 



6 Deriving the PDE 

We wish to derive (1.6) from the identity (2.1). There is a technical issue we need to settle first: in 
(2.2), the function j(x, m, y, n, t) does not have a compact support with respect to (m, n), even if J 
is of compact support. Recall that in Theorem 2.1 we have assumed that J is of compact support. 
Lemma 6.1 settles this issue. 

Lemma 6.1 There exists a constant C4 independent of e such that 

j-T 

(6.1) Eat / e^^'^"^) V a{mi{t),mj{t))V^{xi{t) - Xj{t))mi{t)mj{t)dt < C^. 

Jo 



i,j&I, 



1 



Moreover, 
(6.2) 

rT 

lim supE/vT / e^^'^"^) V a{mi{t),mj{t))Ve{xi{t) - Xj{t))l{mm {mi{t),mj{t)} < L''^)dt = 0. 

Proof. Let us take a smooth function J : M'^ ^ [0, 00) and set 

(6.3) H{x) = co{d) [ Y^^^r^dy 

J \x- y\'^ ^ 

with co((i) = {d — 2)^^uj^^ with ujd denoting the surface area of the unit sphere in W^. Note that 
H >0 and —AH = J. Let ip : (0, 00) x (0, 00) — > [0, 00) be a continous symmetric function and set 

(6.4) X^(q) = £2('i-2) ^ H{xi - Xj)i:{mi,mj) 
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We have 



-Ejv / AcXN{ci{s))ds-EN [ AoXNiq{s))ds = E^Xiv(q(0)) - EivXjv(q(r)) 

Jo JO 

(6.5) < EnXn{cl{0)), 

where 

^o-'^^Ar(q) = -e^^''"^) ^ J{xi - Xj)'^{mi,mj){d{mi) + d{mj)), 

and ^c^Jv(q) = >l(q) + >2(q), with 

^i(q) = -£^('^-^) a(m.^, mj)Ve{xi - Xj)ip{mi,mj)H{xi - a:^) 

l2(q) = e^^''"^^ ^ a{mi,rnj)Ve{xi - Xj)r{xi,Xj,Xk,rni,mj,mk), 



where 



r 5 Xj ^ Xj^^ TTli 5 TTlj , /TZ/i; ) 



rrii + mj 

m 



^^{mi + mj,mk) - i){mi,mk) 



+ 
+ 
+ 



rrii + mj 
rrii + 'Tij 



'^-^'4){mi + mj,mk) - ip{mj,mk) 
ip{mk,mi + mj) - ip{mk,mi) 



ip{mk,mi + mj) - ip{mk,mj) 



H{xi - Xk) 
H{xj - Xk) 
H{xk - Xi) 
H{xk - Xj). 



We consider two examples for ip. As the first example, we choose i^im^n) = mn. This yields 
Y2 = 0. We find that 



(6.6) 



supEjv / Yi{ci{s))ds < E;v^iv(q(0)). 

N Jo 



The hope is that a suitable choice of J would yield the desired assertion (6.1). For this, we simply 
choose J{x) = e~'^A (|) where A is a smooth non-negative function of compact support. We then 
have that H{x) = £^''^3 (|) where AS = -A. As a result. 



(6.7) 
with 



Fi(q)=£''-2 



mimja{mi,mj) 



B{x) = co{d) J 



My) 



|d-2 



dy. 



\x - y\ 

Recall that the support of V is contained in the set y with |y| < Cq. If we choose A so that 

H\y\ < 3Co) < A{y) < l{\y\ < 4Co), 
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then, for |x| < Co, 

B{x) > co{d) [ ^ co{d)C^-'' [ dy =: tq > 0. 

^3Co>|y|>2Co F - y\ ^3Co>|y|>2Co 

On the other hand, if |x| < 5Co, then 

(6.8) B{x) < co{d) [ = lco{d)u;di9Cof. 



'\x-y\<9Co \^-y 

and if |x| > 5Co, then 

B{x) < co{d) 



Ax 



2-d 

J I |2-d 

ay = ci \x\ 

Co>\y\ 

From this, (6.8) and the third part of Hypothesis 1.2, wc learn that the right-hand side of (6.6) is 
uniformly bounded in e. This completes the proof of (6.1). 

As for (6.2), we choose ijj{m, n) = l(m < 6) + l(n < 6). This time we have that Y2 < 0. Such 
a function ^ is not continuous. But by a simple approximation procedure we can readily see that 
(6.5) is valid for such a choice. By the third part of Hypothesis 1.2 on the initial data, we know 
that 

/ / hn{x)h{y)\x — y\^~'^dxdydn < 00. 

From this we learn that 



whence 



lim / / hn{x)h{y)\x — y\^ "^dxdydn = 0, 
5^0 Jo J 



limsupE7vX7v(q(0)) = 0. 



This and (6.5) imply (6.2). □ 

Proof of Theorem 1.1. Step 1. We take a smooth test function J of compact support in 
X (0,00) X [0, cxd) and study the decomposition (2.1). Firstly, we show that the martingle term 
goes to 0. The term Mt is a martingale satisfying 

En[M^]=En {ILY'^ -2Y1LY) {q{t),t)dt = En Ai{ci{t),t)dt + En A2{ci{t),t)dt, 
Jq Jo Jo 

where Ai (q, t) and A2 (q, t) are respectively set equal to 

Ai{q,t) = ^c?(mi)|Ja,(xi,mi,t)|2, 

and 

A2{ci,t) = e^^*^"^^ ^ a{mi,mj)Ve{xi - Xj)J{xi,mi,Xj,mj,tf. 



iei, 



q 
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We can readily show 

(6.9) ^i(q,0 < Cle2('^-2)^d(m,)<C2e'^-^ 

(6.10) Ejv / A2{(i{t),t)dt < csEn / e^^'^'^) V a{mi,mj)Ve{Xi - Xj)dt < c^e"^'^, 

where we have Lemma 4.1 in the last inequality. Prom these inequalities, we deduce that the 
martingale tends to zero, in the e | limit. 

Step 2. We rewrite the terms of (2.1) in terms of the empirical measures. We have that 

(6.11) Y{q{t),t)= / / J{x,n,t)g{dx,dn,t), 

Jo Jilt''- 

and that 

(6.12) ^ (^■^^+Ao^Y{ci{t),t)dt = ^ J^^(^-^^+d{n)A^^ J{x,n,t)g{dx,dn,t). 
Furthermore, by Theorem 2.1 and Lemma 6.1, 

(6.13) / AcY{ci{t),t)dt= f Ti{ci{t),t)dt + Err^{e,L)+Err^{e,S,L), 
Jo Jo 

where T is large enough so that J(-, •, i) = for t >T, the expression r|^(q, t) is given by 

/ a(m,n)U^^{wi — W2)f^{wi,dm;q)f^{w2,dn;q)J{wi,rn,W2,n,t)dwidw2, 
and 

lim sup EAr|£;rr^(e,L)| =0, lim limsupEjv|£;rr2(e, (5, L)| =0. 

L->oo £ S-*0 £_»o 

We note that if we replace f^{w2,dn;ci)J{wi,m,W2,n,t) with f^{wi,dn;q)J{wi,m,wi,n,t), then 
we produce an error which is of order 0{LS~^~^£), which goes to because we send e ^ first. 
As a result, (6.13) equals 

/ / / / P{m,n){g *x ^^)ix,t,dm){g *x ^^)ix,t,dn)J{x,m,n,t)dxdt 
Jo Jr'^ Jl-^ Jl-^ 

+ Err^{£,L) + Err^{e,d, L), 




where 



\imlimsupKN\Err'^{e,5, L)\ = 0. 

5^0 e_>o 



By passing to the limit in low e, we find that any weak limit V is concentrated on the space of 
measures g{dx, dn, t)dt such that. 



hn{x)J{x,n,0)dxdn+ [ f f g{dx,dn,t) ( — + d{n)Ax] J{x,n,t)dt 
Jo Jo Jr^ \ot ) 

(6.14) + / / / / I3{m,n){g *x ^^){x,t,drn){g *x ^^){x,t,dn)J{x,rn,n,t)dxdt 

Jo Vm<* 



+Err'^{L)+Err^{S) =0, 
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where the 7^-expectation of \Err^{S) \ goes to zero as (5 | 0, and the 'P-expectation of \Err'^{L) \ goes 
to zero as L — > oo. From Theorem 5.1 we know that g{dx,dn,t) = f{x,t,dn)dx, T'-almost surely 
and that by (5.2) we can replace g *x ^ with /. Hence 



, hr. 



oo rL 



in{x)J{x,n,Q)dxdn+ I dt f{x,t,dn){— + d{n)Ax)J{x,n,t) 

Jo Jo JRd- \Ot J 

(6.15) +11 I 1 (3{m,n)f{x,t,dm)f{x,t,dn)J{x,m,n,t)dxdt + Err'^{L) = 0. 

Jo Jl-^ Jl-^ JRrf 

It remains to replace and L with and oo respectively. For this, recall that by assumption, 
there exists £ such that J{x, m, i) = if m ^ {£~^,£). Hence, when J{x, m, n, t) / 0, we must have 
that m + n> l"^ and min{m,n} < £. By the first remark we made after the statement of Theorem 
1.1, we know that (3 < a. Prom the second part of Hypothesis 1.1 we deduce that there exists a 
constant C5 = C5{£) such that f3{m,n) < a{m,n) < C572(m)72(n) provided that m + n> £''^ and 
min{m, n} < £. (Here we are using the fact that d{m)'^^'^ (j)'^^'^ is uniformly positive and bounded 
over the interval [£~^/2,£].) On the other hand, we know by part 1 of Theorem 5.1, 

pT p poo /"OO 

/ / j j j2{'n)^2{'m)f {x,t,dm)f{x,t, dn)dxdt < 00, 
Jo J\x\<L-i Jo Jo 

"P-almost surely, where Li is chosen so that the set {|x| < Li} contains the support of J in the 
spatial variable. Prom this we deduce 

pT /• /"OO /"OO 

lim / / / / (5{m,n)f{x,t,dm)f{x,t,dn) 
L-^ooJo J Jo Jo 

1 (max{m, n} > L or min{m, n} < L~^^ J{x,m,n,t)dxdt = 0. 

This allows us to replace and L with and 00 respectively in (6.15), concluding that f{x, t, dn) 
solves (1.1) weakly in the sense of (1.6). □ 

As we stated in Section 1, the family is defined on a compact metric space X which consists 
of measures /x(dx, dn, dt) which are absolutely continuous with respect to the time variable. This 
can be proved by standard arguments. 

Lemma 6.2 Every measure ^ X is of the form fj,{dx, dn, dt) = g{dx, dn, t)dt. 

Proof. Let J^. : M'^ x [0, cx)) M, G N be a sequence of linearly independent continuous 
functions of compact support such that Ji = 1 and the linear span Y of this sequence is dense in 
the space of continuous functions of compact support. Given /x G -Y, it is not hard to show that for 
each k, there exists a measurable function Gj^ : [0, T] ^ R such that \\Gjf, ||loo < Zsup^ ,^ \ Jk{x, n)\, 
and 

/■ /-oo 

Jk{x,n)fj,{dx,dn,dt) = Gjj^{t)dt. 



'ffirf JO 

We wish to define Gj for every continuous J of compact support. Note that each Gj^ is defined 
almost everywhere in the interval [0, 00). For our purposes, we need to construct Gj in such a way 
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that for almost all t, the operator J Gj{t) is linear. For this, let us set Gj = riG jj + • • • + riGj^ 
when J = riJi + ■ ■ ■ + riJi with ri, . . . ,r; rational. The set of such J is denoted by Y'. Since 
Y' is countable, There exists a set A C [0, oo) of Lcbesguc measure, such that for t ^ A, the 
operator J Gj{t) from Y' to M is linear over rationals. By denseness of rationals, we can extend 
J ^ Gj{t) for J G y and t ^ A. For such (J, t), 

/ / J{x,n)iJ,{dx,dn,dt) = Gj{t)dt. 

jR'i Jo 

We then take a point in [0, oo) — A and use Riesz Representation Theorem to find a measure 
g{dx, dn, t) such that 

Gj{t) = / / J{x,n)g{dx,dn,t), 

JlStd. Jq 

for every J eY. Hence 

/ / J{x,n)iJ,{dx,dn,dt) = / / J{x,n)g{dx,dn,t)dt. 
JRd Jo Jw^ Jo 

for every J eY. This completes the proof. □ 



7 Entropy 

In this section, we establish entropy-like inequalities to show that the macroscopic density g is 
absolutely continuous with respect to Lebesgue measure. 

Proof of Theorem 1.2. 

Stepl. Recall that initially we have M particles. We choose /q(o) = {l)---)-^}; a-^d label the 
initial particles as (xi, mi), . . . , {x_\f, mj\f). If a coagulation occurs at time t, one of the coagulating 
particles disappears from the system, and /q C {l, . . . ,J\f} satisfies |-fq(t+)| = |-^q(t)| ~ 1- We write 
A/'(q) = l-fql for the number of particles of the configuration q. Note that A/'(q) takes values in the 
set {1, . . . ,f\f}- We write -F(q, i)^'Ar(dq) for the law of q(t), and define 

HNit) = J F(q,t)logF(q,t) i.^(dq). 

By standard arguments, 

(7.1) ^(*) = / {U^ogF){q,t))F{ci,t)Mdci) = n, + n2, 



where 



= y (Ao(logF)(q,i))F(q,i)i^Jv(dq), 
= J (Ac(logF)(q,t))F(q,t)i/iv(rfq). 
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We have 



iei,_ 

i€lq iG/q 



.F ■ Xi dfN 



X d{mi) — — duN - / X d{mi){d - \xi\'^)F dvjq 
< D \xi\'^F duN, 

•' i€lq 

where we integrated by parts for the second and third equahty, and D is an upper bound for the 
function d(-). To bound the right-hand side, we use the Markov property of the process q(t) to 
write 

Eiv X < Yl l-Ti(0)P + 2d/" Ejv X d{mi{s))ds 



where, in the first inequahty, we used that the coagulation is non-positive, which follows from our 
assumption that a particle, newly born in a coagulation event, is placed in the location of one of 
the departing particles. The second inequality is due to our assumption that Z) is a uniform upper 
bound on d : (0, oo) (0, oo) and to the hypothesis we make on the initial condition. We learn 
that 



(7.2) 



rii < ci(i + l)e 



2-d 



We now concentrate on the contribution coming from coagulations, namely the expression O2. 
This expression equals 



log 



F{Sl.a,t) 



log 



rrii + rrij F(q, t) rrii + mj F(q, t) 



F(q,t) VNidci) 



rrii + "T" 



■F(5i,q,t) + ^F(5?.q,t) 



rui + rrij 



VNidci) 



J X Ve{xi- Xj)a{mi,mj) 

= JY ^e{xi- Xj)a{mi,mj)F{SljCi,t) UN^dci), 

where we used the elementary inequality log a; < x for the second line. To bound this, we first 
observe 

j V,{xi - x,)(27r)-'^/2 (^-^^ dxi < (27r)-'^/2 j V,{xi - Xj)dxi < Ce''-^ 
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We then make a change of variables rrij + rrij i-^ mj. As a result, is bounded above by 

e'^-^y ^pim,)Fici,t)di^Nidq), 



where the function p is defined (II. Sp . 

From the second part of Hypothesis 1.3, we deduce that is bounded by a constant multiple 
of e*^"^. This, the first part of Hypothesis 1.3, and (7.2) yield 

(7.3) HNit) <C2it + l)e'^-^. 

Step 2. Note that by Sanov's theorem, the empirical measure s'^"^ X^j "^(xi.mi) satisfies a large 
deviation principle with respect to the measure vn as £ ^ 0. The large deviation rate function 
I{g) = 00 unless g{dx,dn) = f{x,n)r{x,n)dxdn and if such a function / exists, then 



I{g) = j {flogf - f + l)r dxdn. 



By an argument similar to the proof of Lemma 6.3 of [3j, we can use (7.3) to deduce that if V is 
any limit point of the sequence , then 

j I{g{-,t)) V{dg) <oo, 

for every t. This completes the proof of Theorem 1.2. □ 

8 Appendix: Scaling of the continuous Smoluchowski equation 

We comment on the scaling satisfied by the system (II. ip . under the assumptions that 

d{n) = n-^ 

and 

(8.4) (3{n, m) =TiP + nf, 
with (p,r] G [0,00). Rescaling the equations, 

(8.5) (7n(x,t) = A"/„A7(A^x,At), 
we note that gn satisfies (jl.ip provided that 

(8.6) 1-7(/»-2t = 
and 

(8.7) -a + 7(1 + 7/) +1 = 0, 
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()8.6p ensuring that the free motion term is preserved, (18. 7|) the interaction term. The mass 

poo p 

^/(^) ~ / / fn{x,t)dx dn, 
Jo Jr'' 

which, formally at least, is conserved in time, is mapped by the rescaling to 

(8.8) hg{t) = A"-"'^-2^/i/(At). 
The mass, then, is conserved by the rescaling provided that 

(8.9) a-rd-27 = 0. 

In the critical case, where each of (18. 6p . (18. 7p and ()8.9p is satisfied, we have that 

d/2 - 1 



7 



a 



r/ + (/)d/2 - 1' 
d/2{(l) + ri + l) -2 



■q + cf)d/2 - 1 
and 



i.W) 



2{r] + <l)d/2 - 1) 



In the case that the dimension d = 2, the values 7 = 0, a = 1 and r = 1/2 are adopted, whatever 
the values taken for the input parameters (p and rj. The only critical scaling, then, leaves the mass 
unchanged and performs a diffusive rescaling of space-time. 

Regarding the critical scaling, we recall from Remark 1.2 of [6] that the condition rj + (p = 1, 
which is a natural transition for the rescaling gn (as is apparent from (j8.10p ). represents the limit 
of the parameter range for which uniqueness and mass-conservation of the solution of p.ip are 
proved: indeed, the condition required by (Bj is ry-|-0 < 1, along with some hypothesis on the initial 
data. 

Do we expect the complementary condition rj+(j) > 1 to have physical meaning? To consider this 
question, we take positive and fixed (p and rj, and consider the rescaling (j8.5p under the constraints 
(j8.6p and (|8.7p . Seeking to understand the formation of massive particles, rather than spatial 
blow-up, we fix r = 0. We are led to 

(8.11) 7 = r' 
and 

(8.12) a = l + ^-^. 



Returning to (j8.5p . a self-similar blow-up profile is consistent with the scaling 
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given by A = provided that its mass (18. 8p does not grow to infinity as A — > 0. We have set 
r = 0: as such, the condition that ensures this is a — 27 > 0, which, by (18. lip and (I8.12p . amounts 
to the inequahty (/> + ??> 1. 

We conclude that considerations of scahng would in principle permit a blow-up in the equations 
in the mass variable under the condition that ij + (p > 1. The blow-up we considered is in a low 
A limit, which corresponds to heavy mass at late times: as such, it should be considered not as 
a gelation, in which particles of infinite mass develop in finite time, but rather as the appearance 
of populations of arbitrarily heavy particles at correspondingly high time-scales. Expressed more 
precisely, the weak form of blow-up considered is the statement that, for each K strictly less than 
the total initial mass f^a mfmix, 0)dxdm and any mg G M"*", there exists t G [0, 00), 



.13) / / mfm{x,t)dxdm > K. 

Jmo Jr'' 



(This condition is correct in the absence of gelation. Gelation would remove mass from all finite 
levels. Note also that the absence of fragmentation in (jl.ip means that, in fact, (I8.13P implies the 
stronger statement that most of the mass accumulates in arbitrarily high levels at all sufficiently 
late times.) In dimension d > 3, (1.11) of Theorem 1.1 in [6] shows that the discrete analogue of 
([8^3]) fails iir] + (j)<l. 

A parallel may be drawn between the Smoluchowski PDE and the non-linear Schrodinger equa- 
tion. Consider, for example, a solution of cubic defocussing NLS, ti : M-^ X R+ ^ C of 

d 

(8.14) i—u-Au = -\u\'^u, 
^ ' dt ' ' 

may be written in Fourier space as 

(8.15) i—u—\^\^u = —j j — ri)u{a)u{ri — a)drjda. 

We see that the mass variable in (jl.ip may be viewed as analogous to the frequency variable in 
(j8.15p : the non- linear interaction term in each case is a type of convolution. Pursuing the analogy, 
the quantity |||Vn||2 + jH^iHl is formally conserved in NLS, as is the mass Jg°° jj^^j mfmdxdm for 
the Smoluchowski PDE. For NLS, the term weak turbulence refers to the growth to infinity in time 



for some s > 1, a circumstance that is anticipated in ()8.14p in a periodic domain. (See Section II. 2 
of [1] for a discussion.) The counterpart of weak turbulence for the system (jl.ip is 

I I rrf fm{x,t)dxdm —>■ 00 as t —>■ 00, 
Jo JR'i 

for some r > 1. (Note that (|8.13p implies this statement for every r > 1 on a subsequence of times.) 
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Comparing the system (ll.lh to its spatially homogeneous counterpart, given in the discrete case 
by {on : [0, oo) — > [0, oo) : n G N} satisfying 

d oo 
(8.16) ^ ^ ^ ~ '^)"™(*)'^«-"^(*) ~ 2 ^ n)am{t)anit), 

m=l m=l 

we see the stabilizing role of diffusion: for example, it is easy to see that, taking /3(n, m) identically 
equal to a constant in ()8.16p ensures the analogue of ()8.13p . while we have seen in the spatial case 
that scaling arguments do not disallow ()8.13p under the condition that + 1. 

Regarding the prospect of proving mass-conservation for at least some part of the parameter 
space where (j) -\- r] > 1, we comment that, in [B], hypotheses of the form P{n,m) < n'^ + 
were used. It may be that, if l3{n, m) < n^^'^ + m^~^^ or /3(n, m) < fi'^/'^+<^jji^/'^+'^ (with e > a 
small constant), but /3 is permitted to have space-time dependence subject to such a bound, then 
gelation is more liable to occur. As such, an argument for mass-conservation would have to exploit 
the assumption that P{n, m) is constant in space-time, in a way that those in [6] did not. 
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